Nonlinearities of the Vlasov equation in the spinodal region 
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Abstract 

We estimate the importance of the nonlinear terms in the Vlasov equation for the development of 
the unstable modes. The results allow to identify the region of wavelength where the linear evolution 
is justified. 
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Recently the Vlasov equation has been applied for the description of the dynamics of the fragment 
^ | formation in the nuclear matter unstable against the spinodal decomposition |l|, |^, In the linear 

approximation a perturbation of the translationally invariant distribution is unstable. Unstable modes 
will develop exponentially in time leading to the formation of density inhomogeneities and fragmentation 
pi Q. Since the growth rate depends on the wave- vector, the wave- vector with the largest growth rate 
T/t will dominate at large times. One expects the formation of fragments of the size corresponding to 
the most unstable modes 1/k Q. If the linear regime can be trusted in the development of the spinodal 
decomposition then the fragment size distribution should show an excess of fragments of intermediate 
sizes. The fragment size distribution would be different if at the fragmentation time many different modes 
are excited due to the nonlinear evolution of the Vlasov equation. Of course, in realistic simulations initial 
fluctuations and a noise term must be taken into account. This determines the strength by which different 
modes of the linearized Vlasov equation would be excited. 

The nonlinearities in the Vlasov equation are non-negligible at large times as the deviation of the 
unstable modes from equilibrium becomes important. The existing estimates of this effect are based on 
the numerical results for the 2-Dimensional Vlasov equation [|[ In this letter we shall estimate the 
importance of the nonlinear effects in the Vlasov equation analytically. We shall study the case where 
at initial time t = one or two unstable modes are perturbed. One unstable mode engenders nonlinear 
effects at double wave-vector. Two unstable modes can lead to nonlinearities also for small wave-vector. 
The estimate will allow to define the regime of wave- vectors where the linear approximation can be valid. 
We start with the usual Vlasov equation : 

m £ ,P) +W x f(t,x,p) - V x U(x)V p .f(t,x,p) = , (1) 
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where U{x) is the density dependent mean field. In the following we take a Skyrme type mean field 
potential folded with a Gaussian : 

U {x )=l^y 9{x -y)(A^ + B(^-Y). (2) 
J v Po v Po ' ' 

Expanding to second order in Sf(t,x,p), the deviation from the homogeneous solution, and Fourier 
transforming in x, 8f(t,k,p) — J d 3 x 5f(t 7 x,p)e~ lxk , we obtain : 



2 Sp 2 de J (2tt) 3 
i f dH 6U k -i 



2 J (2tt) 3 Sp 



8p(t,k-l)QL-l)V p 6f(t,k,p) = 0, (3) 



where hq is the equilibrium momentum distribution. The first three terms of the above equation are the 
linear approximation to the Vlasov equation. The remaining two terms are the second order expansion 
in Sf and represent a mode-mode coupling for the modes of the linear equation. The actual solution of 
the second order equation is as difficult as for the full Vlasov equation. However, to make an estimate 
of the nonlinear effects it is enough to take for 6f(t,k,p) and Sp(t,k) in the nonlinear terms only the 
contribution from the unstable modes to the evolution of an initial phase space density perturbation at 

Sf + (t,k,p)= S ^ d ^-^5 P +(k)e^, (4) 
dp oe kv — il i 

and 

Sp+(t, k) = 6p+(k)Q(k max - \k\)e Tkt , (5) 
with k max being the maximal wave-vector with imaginary frequency 0] and 

jr + m f d 3 p 6f(0,k,p) 

* pW = iWM' (6) 

Inserting the forms (|B|) and (^|) in the mode-mode coupling term in eq. ([}]) and taking the one-sided 
Fourier transform : 

/>OG 

5f(u,k,p)= / dtSf(t,k,p)e iut , (7) 
Jo 

with Im uj > 2T rnaxl one obtains : 

i(kv -u)Sf(u), k,p) - i^^kv^-5p(u,k) 

op Oe 

i S 2 Uk. dn f d 3 l 
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iuj + Tk-i + Ti 

i f dH SU k -iSUi s ,„ nr 1 



-Sp + {k-l)Sp+{l)- 



2j (2tt)3 Sp Sp p v /p w z^ + r fc _, + r ; 



Av-iTi de J 

Dividing by kv — u and integrating over p an equation for the density perturbation is obtained : 

-iG(u),k) 
Sp{uj,k) = — - — — 
e{u 7 k) 
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1 S 2 U k f d 3 p kv dn 



d 3 l 

2e(uj,k) Sp 2 J (27r) 3 kv-w de J (2rf 



■Sp + {k-l)5p + (l) 



1 



iuj + r fe _; + Ti 
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2e(w,fc) 
d 3 p 



(2tt) 3 kv - uj 



dH SU k ^ SU, 
(2tt) 3 Sp Sp 

' (k - 1)V P 



S P + {k-l)Sp + {l) 



1 



lv 



dn \ 



lv - de / 



with 



G(uj,k)= J 



and 



e(oj,k) = 1 



9C4 
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d 3 p 6f(0,k,p) 
(2tt) 3 kv-w 

d 3 p kv 



dn 



(27r) 3 kv-w <9e ' 

The time dependence of the density perturbation Sp can be found using the inverse transform 



Sp(t,k) 



^Sp(uj,k) 



the integration path in the inverse Fourier transform laying above any singularity of the integrand. 
As a first case we take only one unstable mode : 

p+(k) = (2Tr) 3 S 3 (k-k )SA koP . 

The nonlinear term is nonzero only for the mode with the doubled wave-vector 2fcp : 



0) 
(10) 
(11) 

(12) 
(13) 



Sp(u>, k) 



-(27r) 3 S 3 (k-2k )5Aiy 



e(u),k) 



1 



feu 



i s 2 u 2kl 

8tt 2 Sp 2 

su k0 



k 2 
ft-o 

~4tt 2 



Sp 
2T k 



k v{uj - 2iT ko ) 2 
k v(u) - 2iT ko ) 




(14) 



The integral (|T^) can be calculated closing the integration path in the lower half-plane. Its value is 
determined by the singularities of the integrand. The inverse susceptibility l/e(u>,k) has two poles 
corresponding the solutions of the dispersion relation : 

e(uj,k)\ k = 0. (15) 

The integrand also has a pole at uj = 2iT ko and a cut on the real axis in uj [^|. The contribution from 
the pole at 2iT ko is dominant at large times. Even if the the solution of (^|) is imaginary we always have 
2T ko > T2 ko ■ The cut contribution is always bounded and can be neglected at large times in comparison 
to the growing components. The result is : 



Sp(t, k) 



(27r) 3 S 3 (k-2k )SA 2 p 2 



i(2iT ko ,2k ) 



k 2 



8tt 2 V Sp 



SU, 



feo 



p 2 dp 



5 2 U 2ko , 25U ko 
Sp 2 Sp 

rfe dn 



((M 



r 2 fc0 ) 2 de 



(2n) 3 S 3 (k-2k )F 2ko (t) 



(16) 
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In Fig.l is shown the ratio of the amplitude F2k a (T) of the mode at 2fco to p at the instability time 
T which is defined as \5Ak exp(Tk T)\ ~ 1 jj], Q. The result shows that the long wavelength modes 
are always nonlinear. Before the time the instability drives the mode fco to large values the mode with 
double wave-vector becomes larger. However, the amplitude of the linear modes with large wave-vectors 
(k > Afm^ 1 ) is large (\SA k exp(TfcT)| ~ p) before the nonlinearly driven mode at 2fco becomes important. 
In particular in our case the most unstable mode is in the region of wave- vectors, where the linear regime 
is valid up to the instability time T. Similar observation were made in the 2-Dimensional numerical 
solutions of the Vlasov equation, the small wave-vector modes become nonlinear at smaller amplitude 
than the large wave- vector modes ||. In Fig. 2 is shown the growth time for the amplitude of the mode 
to reach p. We see that for the long-wavelength the growth time for the nonlinear mode at 2fco is smaller 
than the growth time for the linear mode at fco. The unphysically long values of the growth times depends 
of course on the initial perturbation chosen (here 8 A — 1/50). 

The nonlinearities become extremely important for long wave-length (Fig.l). However in this limit 
one should take into account both the pole at 2iT ko and at i^2k , since the two poles merge in that limit 
(T k ~ k). The time dependence is now different : 



(27r) 3 (5 3 (fc-2fco)M2 p 2 9P , 
S P (t,k) = [ L- \ °> koP te^ot 
le (w,2fc )^ =2 ;r fco 



S 2 U 2ko , 25U ko 
Sp 2 Sp 



' 1 ' p dp-. 



8tt 2 V Sp J J y ((k v) 2 +T 2 o ) 2 de 
= (2n) 3 6 3 (k - 2k )L koP SA 2 ko t e 2Ft . (17) 

Now, the amplitude F2 ko (T) at the instability time T depends on the value of the initial perturbation 
SA ko : 

F 2k0 (T) = -^H\SA k0 \)p. (18) 

1 k 

This shows that it is not the same to use larger initial fluctuations evolved on shorter time to using small 
initial fluctuations evolved for longer times. Although, the linear evolution would give the same result 
|(SAexp(rT)| = 1, the nonlinear modes at double wave- vector would have a different strength. In long 
wavelength limit the procedure Q of putting instead of the noise term in the Vlasov equation a stronger 
initial perturbation gives different results due to nonlinearities of the Vlasov equation. However, the 
most unstable modes are behaving according to eq. ([l6|). Moreover, the numerical coefficient is such 
that the nonlinearities have no time to built up until the amplitude of the most unstable modes becomes 
of the order p. The results from the correct limiting expression for the limit T(k) ~ k (eq. |l^) is also 
shown in Fig. 2. The correct dependence of the growth time T growt h ~ 1/fc in the long wavelength limit 
is recovered. 

Another case is when two unstable modes develop. We shall study the case when the wave- vectors of 
the unstable modes are in opposite direction : 

p+(k) = (2tt) 3 (S 3 (k - l)5A l + 5 3 {k + k ~ l)SA ko -i) p , (19) 

with kg < I and / corresponding to a strongly growing mode. Besides the contributions at doubled 
wave-vector 21 and 2(fco ~ we find a nonlinear mode located at fco. The mechanism responsible for the 
growth of the mode at fco is different from the simple wave- vector doubling. As previously the dominant 
mode is given by the pole of @ at i(Ti + Tk -i), leading to the behavior : 

F fco (t)~M fcD _,^c( r '+ r * D -0* . (20) 

We do not quote explicitly the lengthy formula for the coefficient, given by residue at the pole of (^|). In 
Fig. 3 is shown the ratio \F ko (T) / F ko _i(T)\ at the time T when the amplitude of the most unstable mode 
I ~ .7/m _1 reaches p. For |fco — l\ < .3/m _1 the linear mode at fco — I is weaker than the nonlinear 
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mode at ko > .3/m _1 . Similar nonlinear effect were observed in numerical simulations 0. In the case of 
a strongly excited mode at relatively small wave-vector and a weakly excited mode with strong growth 
rate several modes at intermediate wave- vectors appear, leading to strong nonlinearities. Of course if the 
most unstable mode I is strongly excited it will dominate both the linear mode at ko — I and the nonlinear 
mode at ko- When ko — > I the contribution from the pole of (Q) at Tk cannot be neglected. This changes 
the behavior of the time dependence of the nonlinear mode from (|2(]) to a form analogous to ( p7| ) . Thus 
the singularity in Fig. 3 when ko —>■ I = .7 fm^ 1 is spurious. 

In summary, we have studied the importance of the nonlinearities in the Vlasov equation for the 
development of the spinodal instabilities. The mechanism of wave-vector doubling is important only in 
the regime of wave- vector k < .35/m -1 (this value is determined mainly by the range of the interaction). 
For the most unstable modes the nonlinear effects at double wave-vector are small. For the case of 
two exponentially growing modes with wave-vectors in opposite direction, we have found that nonlinear 
modes at intermediate wave-vectors can be excited. Thus, a small excitation of a rapidly growing mode 
~ HjmT 1 together with a strong excitation of the mode ^.3/m _1 leads to strong nonlinear effects at 
wave- vectors ~ .4/m _1 . If for a given wave- vector the nonlinearities are found to be important, then the 
perturbative approach leading from (^|) to (|§|) becomes invalid at some time smaller than the instability 
time. From the one and two mode cases studied in this letter any initial one-dimensional perturbation 
of the homogeneous solution can be constructed. In the second order in 5f(t,k,p) the corresponding 
nonlinearities from mode-mode coupling can be calculated. Those estimates explain the observations 
made in numerical simulations of the 2-Dimensional Vlasov equation 0, || and are the first estimate of 
nonlinear effects in 3-Dimcnsions. 
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Fig.l 

The ratio of the amplitude \F2k(T) \ of the nonlinear mode at double wave- vector 2k to the average density 
p as a function of the wave-vector at the instability time T = ln(l/|(L4fc|)/IV The calculation is done 
at zero temperature with the parameters of the interaction A — ~356MeV, B = 303MeV, a = 7/6, 
p = pa/3 and the range of the Gaussian equal to .9/m. 
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Fig.2 

The growth time T grow th for the amplitude of a mode to reach the average density p with 5Ak = 1/50, 
as a function of the wave- vector. The solid line is for the nonlinear mode at 2k, the dashcd-dotted line is 
for the linear mode at k and the dotted line is for the nonlinear mode at 2k in the regime Tk ~ k. The 
parameters used are the same as in Fig.l. 
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Fig.3 

The ratio of the amplitude Fk(T) of the nonlinear mode at k to the amplitude Fk-i(T) of the linear mode 
at k — I as a function of the wave-vector at the time T = ln(l/\SAi\)/Ti (I = .7/m -1 ), i.e. when the 
amplitude of the other linear mode reaches the average density. The parameters used are the same as in 
Fig.l. 
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